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Within Ginzburg-Landau theory, we study the structure of a magnetic vortex in color-flavor locked 
quark matter. This vortex is characterized by winding of the 517(3) phase in color-flavor space, as 
well as by the presence of a color-flavor unlocked condensate in the core. We estimate the upper and 
lower critical fields and the critical Ginzburg-Landau parameter that distinguishes between type I 
and type II superconductors. 
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I. INTRODUCTION 



o : 
o . 

The possibility that dense quark matter occurs in neutron stars has been considered for the past three decades. The 
possible presence of a diquark condensate in such quark matter Q has generated new interest in such matter. Among 
various condensates, a color-flavor locked (CFL) state Q, in which all three color and flavors are gapped, has been 
studied most intensively, along with a two-flavor pairingstate 0- In the weak coupling (high density) and massless 
I ■ limit the CFL state is most stable at zero temperature Q and near the transition temperature a feature ensured 
\ by the quark-quark attraction in the color antitriplct channel due to one-gluon exchange. Although the massless, 
weak coupling results for diquark condensation are modified by effects of strong coupling and nonzero quark masses 
in a regime as might be encountered in possible quark matter in neutron star interiors |fj, the CFL phase has to be 
^ . reckoned with as a reference state in such a regime. 

Response of the CFL condensate to ma gne tic fields, as would be experienced by quark matter if present in neutron 
stars, was considered in Refs. IMHISlEllfl- Because of nonvanishing color and electric charge carried by diquark 
■ pairs, electromagnetic superconductivity and color superconductivity occur in the CFL phase. These two phenomena 
t—{ [ are not independent in the sense that transverse photon and gluon fields are mixed with each other. One of the 
• resultant mixed fields is freely propagating, while the other is Meissner screened and thus massive. The massive 
mixed field and the winding of the SU(3) phase in color-flavor space are essential to the structure of a magnetic 
vortex in the CFL phase. This point was not taken into account in Ref. The criterion of whether or not the CFL 
O |' state can allow magnetic vortices to form (type II or type I) has been described by Giannakis and Ren [l(| in terms 
of the energy needed to form a planar surface separating the normal and superconducting material, but the detailed 
structure of the magnetic vortex has yet to be determined. 

In this paper we examine the properties of magnetic vortices in the CFL phase within the framework of Ginzburg- 
Landau theory. These vortices are characterized by three different fields: the CFL diquark field, the massive photon- 
gluon mixed field, and a color-flavor unlocked diquark field coupled to the mixed field [T(il ]. In the London limit 
where the wavelengths of distortions of the condensate are long compared with the size of the vortex core, the 
supercurrents surrounding the core can be written in terms of the color-flavor SU(3) phase of the CFL condensate 
and the mixed field. From such supercurrent structures we derive the lower critical field, H c \. We next examine the 
vortex structure near the core, in which the color-flavor unlocked condensate plays a role, and derive the upper critical 
field, H C 2- Finally we revisit how to distinguish between type I and type II superconductors from the viewpoint of the 
intervortex interactions. Throughout this paper, we consider a system of three-flavor (uds) and three-color (RGB) 
massless quarks at temperature T and baryon chemical potential /i, and use units K = c = 1. We assume that the 
Fermi momentum is common to all colors and flavors since in the massless limit of interest here, it is the case with 
both the normal and CFL phases of equilibrated neutral quark matter. 



II. GINZBURG-LANDAU THEORY 

We first review the Ginzburg-Landau theory relevant for the description of the CFL phase [f|, 0, . While this 
theory is strictly applicable near the critical temperature T c and for inhomogeneities of wavelengths longer than the 
zero temperature coherence length £o ~ T~ , it is nevertheless useful in describing the structure of a possible vortex in 
the presence of external uniform magnetic fields. This theory focuses only on the pairing channel that first prevails in 
a normal quark-gluon plasma as the temperature goes down. In the weak coupling (high density) regime, this pairing 
channel is predicted to have even parity, zero total angular momentum, and same chirality, and to be in a color and 
flavor antitriplet state. As we shall see, the CFL phase belongs to this pairing channel. We assume that this channel 
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has the highest transition temperature also at lower densities. This pairing state is characterized by a complex 3x3 
gap matrix, (d Q )j(r), in color-flavor space [5|, where a (i) is the color (flavor) other than two colors (flavors) involved 
in Cooper pairing. This gap is defined on the mass shell of the quark quasiparticle of momenta on the Fermi surface, 
but remains dependent on the center-of-mass coordinate r of the quark Cooper pair. 

In the presence of an external magnetic field, H cxt = V x A oxt , one can write down the thermodynamic potential 
density difference Afl — f2 s — fl n between the superfluid and normal phases near T c as @ 

AO = aJ2 KI 2 + A£ KI 2 ) 2 \< ■ d b\ 2 + 2R t E \( D M 2 + \ G lm G Tm + -F lmFlm - i|H cxt | 2 . (1) 



a b 



Here we have neglected the response of the normal component to H oxt . The parameters a, and /?2 characterize the 
homogeneous part of the free energy, while Kt is the stiffness parameter, controlling spatial variations of the order 
parameter. Since (d a )i is antisymmetric in color and flavor space, the covariant derivative Di is 

(D,d)„ = d,d a - l -gA?(\ a d) a - ieAi(Qd) a) (2) 

where the A" are the Gell-Mann matrices. Gf m (Fi m ) is the spatial part of the gluon (photon) field-strength tensor, 
Q = diag(2, — 1, — 1)/3 is the electric charge matrix in flavor space, and g is the QCD coupling constant. Here we 
have adopted the convention in which A 8 = diag(2, — 1, — 1)/\/3 and (d a )j is defined in terms of the adjoint spinors, 
rather than the spinors 

The parameters a, fli, 02, and Kt in Eq. (JIJ can in principle be determined as functions of T and /i. In a general 
Ginzburg-Landau framework, however, these parameters are left unknown. The parameters have been determined 
only in the weak coupling limit in which the pairing interaction is controlled by Landau-damped magnetic gluons, 
andareHEl 

a = 4iV(/V3) In [JV (3) 



Tc, 

U = h = *K, =^^j2 N M 3 )> ( 4 ) 



with the zeta function £(3) = 1.202 • • •. 

In homogeneous quark matter in weak coupling, the energy difference Q has a minimal value [f| just below T c for 
a CFL order parameter, i.e., 

{d a )i = U ai K A , K A = e^°\K A \, (6) 

where U = exp(i\ a ip a /2) represents SU(3) rotations in color-flavor space. Generally, in a homogeneous situation, the 
CFL state is the most stable as long as 3/?i + (i 2 > and /3 2 > 5j. In terms of a, f3\, and /3 2 , the minimizing value 
of Ik^I is 

\6/3cfl/ 

where /3 C fl = Pi + 

We can now derive the field equations characterizing the CFL state under constant external magnetic field, H cxt . 
In this situation, as clarified in Ref. [l(j> the electromagnetic field, which is mixed with the gluon field of a = 8, is 
coupled with the color-flavor rotation exp(zA 8 (/?g/2) of the gap matrix (d a )j. Interestingly, the solution to the field 
equations requires a nonvanishing color-flavor unlocked diquark component proportional to [A 8 exp(iA 8 ys/2)] a i, in 
addition to the CFL component proportional to [exp(iA 8 (pg/2)] a i |lOj : addition of any other components would simply 
increase the total free energy of the system. It is thus convenient to introduce the ansatz for the pairing gap [To| : 



(d a )i(r) = 5 aR S iu (j>(r) + (5 aG Sid + S aB 5i S )x(r)/V2 



a iA-„.(r)/ 2l |<Hr)|+V2| X (r)| 



mv)\-V2\ X (v)\ 

71 ' (8) 
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with (j>(v) = exp[i</?8(r)/v / 3]|0(r)| and x(r) = exp[— iyg(r)/2v / 3]|x( r )li an d the photon-gluon mixed fields [T^: 

A _ V3gA~2eA s ^ _ V3gA s + 2eA 

6ffCFL * 6gCFL 

where gcFh = y/3g 2 + 4e 2 /6. Then, we can express the thermodynamic potential density difference as 

AO = a(|0| 2 + |x| a )+)9i(H 2 + |x| 2 ) 2 +A(H 4 + WV2) 

+2X T | (V - 2 lffCFL A 8 ) 0| 2 + 2^ T | (V + igcFLA 8 ) xf + i(|B| 2 + |B 8 | 2 - |H cxt | 2 ), (10) 

where 3 — V x A and 8 8 = Vx ,4 8 . By extremizing this free energy difference with respect to the pairing gaps and 
gauge fields, we obtain the gap equations, 

- 2K T {V ~ 2ig CFh A 8 ) 2 (j ) + a<t> + 2/3 1 (|0| 2 + \ X ?)4> + W2\4>?<t> = 0, (11) 

- 2K T (V + ig C FLA s ) 2 x + a X + 2/3i(|0| 2 + |x| 2 )x + /? 2 |x| 2 X = 0, (12) 

and the classical Maxwell equations for the photon-gluon mixed fields, 

V x (V x A & ) = Re{UK T g CFh [-24>*(V - 2ig CFh A s )0 + x*(V + z.g CF L.4 8 )x]} 

= J\ (13) 

Vx(VxA)=0. (14) 

From Eqs. (|13fl and (|14f> we find that in the CFL state (\<f>\ = \x\/V% = I k a|)> the mixed field A s is Meissner screened 
within a London penetration depth, 

CFL 2^6K T gcFh\KA\ ' 

whereas the mixed field A is freely propagating. 

The thermodynamic critical field, H c , associated with B = V x A is the field at which the Gibbs free energy density 
of the normal state, G n = Q n — H cxt • B, drops to that of the superconducting state, G s = 0= — H cxt • B. Since B s is 
screened out in the CFL phase in bulk, one can obtain the Gibbs free energy difference as 

a 2 p 2 

AG = G s -G„ = --= + —^5 — |H cxt | 2 . (16) 

4/3CFL ioffcFL 



Then, 



He = , (17) 

^ e sCFL^CFL 



2i^ 1/2 



where 

^={laf) (18) 

is the Ginzburg-Landau coherence length. As discussed in Ref. extrapolation of the weak coupling expression 
for H c to low densities and temperatures indicates that the critical field H c is typically ~ 10 19 G, several orders of 
magnitude larger than canonical neutron star surface fields ~ 10 12 G. 

We conclude this section by asking whether or not the CFL condensate allows magnetic vortices associated with 
the field A 8 to form, i.e., whether it is type II or type I. In a standard method for distinguishing between type I and 
type II superconductors, one calculates the energy per unit area, a s , needed to form a planar surface separating the 
normal and superconducting material. At the thermodynamic critical field, (|17fl . applied parallel to the surface, the 
surface is in mechanical equilibrium. For a surface perpendicular to the z axis, the surface energy can be written as 
the integral over z of the difference between the total Gibbs free energy density and the value of \z\ — > oo 0: 

/oo 
dz[An(z)\ Hoxt ^ Hc - H c (\B(z)\ - H c )], (19) 
-oo 
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Because of the presence of the color- flavor unlocked component near the surface, the system remains type II, i.e., 
a s < 0, even when the Ginzburg-Landau parameter, kcfl = Acfl/£cfLj becomes less than l/v2 by a small amount 
[To|. This is a contrast to the case of ordinary superconductors in which the Ginzburg-Landau parameter is l/\/2 
when the interfacial energy vanishes [l4j . Detailed calculations 01 show that for /3i = @2 the system remains type 
II for values of kqfl down to ~ 0.589. As we shall see in Sec. II VI however, the properties of intervortex interactions 
imply a significantly narrower type II region. 



III. MAGNETIC RESPONSE IN THE LONDON LIMIT 

We now examine the response of the CFL condensate to external magnetic fields by deriving the supercurrent S 8 
in the London limit in which a length scale of the spatial variation of the condensate is large compared with the 
coherence length £gl, corresponding to the core size of a magnetic vortex. In this limit, one can ignore the variation 
of the magnitude of the pairing gaps, <p and as compared with the variation of their phases We may thus derive 
J 8 by substituting the ansatz © with \4>(y)\ = |x(r)|/\/2 = into Eq. {T3J). The result reads 

J 8 = 4K T g C FL\KA\ 2 (V3V(p 8 ~ 6 5c flA 8 ). (20) 

We thus find that this supercurrent can be induced by the gradient of the color- flavor SU(3) phase, cp$- This feature 
is essential to magnetic vortices, but was ignored in Ref. 9]. From expression (|20|l we can examine the response 
properties of the CFL condensate. 

We first consider the response to weak uniform magnetic fields. As discussed in Ref. [9j , the response is characterized 
by imperfect diamagnetism, i.e., partial Meissner screening. As long as g > e, most of the field freely propagates in 
the form of B, whereas the rest is included in B 8 and hence screened on a length scale of Acfl- 

For high fields, it is interesting to consider the possible presence of magnetic vortices. We find from Eq. (|20[1 that 
vortices can appear in such a way as to satisfy the flux quantization condition, 

idl-(A 8 + \ 2 CF ^J 8 )^^ = 4> S n 1 (21) 

J ffCFL 

where the integration is performed around a closed loop surrounding the vortex line, </>g is the flux quantum, and n 
is the winding number of the vortex. Whether vortices actually occur in equilibrium depends on the value of kcfl, 
as will be discussed in Sec. IIVI 

In the London limit we can analyze the structure, far away from the core, of a singly quantized (rt = 1) vortex and 
then estimate the lower critical field, H c i, above which a single vortex starts to appear. Let us take the vortex to be 
aligned along the z axis and set <p$ — 2\f3ip, where ip is the azimuthal angle around the line. Then, the field equation 
(|13fl reduces to the London equation, 

1 2n 

# 8 + AcfiV x (V x B 8 ) = —j= V x Vtps = 5(x)S(y)z. (22) 

2v 3<7CFL 3CFL 

Note that this equation leads to Eq. I|21|l with n = 1. 

As in ordinary type II superconductors |14|. one can solve the London equation (|22|l and obtain the supercurrent 
density flowing in the azimuthal direction around the line as 

08 

^— T5 <Pi for £cfl < r < Acfl, 



2ttA cfl V 2r 



CFL 



'•P 



Acfl 2r 



^FL 



-r/A CFLj for r> Acfl- 



The decrease of this supercurrent at a scale of Acfl suggests that the field B s associated with the vortex can penetrate 
only up to such a scale, beyond which A. 8 ~ (l/gcFLr)<p and thus B 8 ~ 0. Combining Eqs. I|21|l and (|23(l . one can 
calculate the vortex line energy per unit length, which is composed of the sum of the magnetic and flow energies, as 

T L = . ff ln^cFL- (24) 

At the lower critical field H c \, the line energy is balanced by the energy gain due to the magnetic induction B 8 = B 8 z. 
The energy gain per unit volume is -Emag = —B 8 T-l 8 , where Tt 8 = (e/3<7cFL)-ffci is the critical field associated with the 
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massive photon-gluon mixed field [see Eq. (JSJ]. Integrating over the volume including the vortex, one can estimate a 
total energy gain per unit length as —<ps( e /^9cFL)H c i. Then, 

3 

H ^ = 7TT2 l nK CFL- (25) 

Note that the logarithmic factor in this expression contains uncertainties of order unity. 



IV. STRUCTURE OF MAGNETIC VORTICES 



As shown in the previous section, the analyses in the London limit are useful in clarifying the global properties of 
the magnetic response in the CFL phase. However, it is important to examine the vortex structure near the core, 
e.g., in evaluating the upper critical field H C 2 where the vortices essentially fuse and the system becomes normal. In 
this section we thus go back to the Ginzburg-Landau equations and examine their solutions. Using these results, we 
revisit the problem of how to distinguish between type I and type II superconductors. 

It is useful to start with the linearized gap equations with respect to 4> and \ : 

- 2K T {V - 2ig CFh A s ) 2 (j) + a<f> = 0, (26) 



- 2K T (V + igcFhA 8 ) 2 X + a X = 0. (27) 

This is because equations (|26|) and 127fl can describe the situation near H C 2 in which the pairing gaps are suppressed 
in magnitude. The critical condition for the existence of nontrivial solutions to at least one of these equations is 
satisfied when 5cfl|# 8 | = \(x\/2Kt- In this situation, a nontrivial \ solution to Eq. (|T7Jl occurs, while = 0. Using 
Eq. © we thus obtain 

H c2 = (28) 

?CFL 

Combining Eqs. IT7jl and jSRjl. we find that H C 2/H C = 2kcfl- When H C 2 > H c , i.e., kcfl > 1/2, magnetic vortices 
are expected to occur. This condition in fact disagrees with the criterion for type II superconductors as derived in 
Ref. M 

from the sign of the normal-super interfacial energy a s (see Sec.|nJ. Note that in ordinary superconductors 
[l4j . the condition H C 2 = H c agrees with the condition a s = 0. The difference between these two cases arises from the 
fact that in the CFL case, magnetic vortices are characterized by the color-flavor locked and unlocked condensates, 
while in the ordinary case, they are characterized by only a single complex scalar condensate. 

This two-field nature of the CFL vortices can be seen by substituting cf>(r) — \ka\ exp(2irup)u(r) , x( r ) = 
v^I^aI exp(— in<p)v(r), and A, s (r) = A^(r)<p into the field equations i|12[l . and JT3J|. Then we obtain 

/ (ft I d \ /n \ 2 

- 2A> ( — + -— J u + 8K T ( - - g CFL Al) u + au + 2/3 1 \k a \ 2 (u 2 + 2v 2 )u + 2/3 2 M V = 0, (29) 



2R t ( TT + -T ) v + 2Kt I" " 9cflA*) v + av + 2(3 1 \k a \ 2 (u 2 + 2v 2 )v + 2P 2 \ka\ 2 v 3 = 0, (30) 
ar z r dr J \r * J 



rAl = , + « 2 ) (- - 9cflK) ■ ( 31 ) 



did 

drrdr'"* v 3g C Fh^cFh 



In the limit r — * 0, the normalized gaps u and v and the gauge field ^4* are suppressed in magnitude. Consequently, 
the gap equations reduce to 

d 2 1 d \ /n\2 , 

— + --U + 4 - u^0, (32) 
dr z r dr J \ r J 

d 2 ^1 d\ + f n \ 2 o (33) 
dr 2 r dr J \r ) 
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From Eqs. ^ and we find that 

u~>dr 2n , v^C 2 r n , (34) 

where C\ and C 2 are real constants that can be determined by connecting the solution near the center (r = 0) with 
that outside the core. We can then calculate the asymptotic behavior of A^ from 

did ^ 8 1 ^„ , . ,,n 

drrdr v 35cflAc FL 

the result reads 



: r < "> o J3— ( 2 « 3 + v2 )r; (35) 



/|8 -^-B s r ^ r 2n + 1 CifiS 

2 12pCFLA CPL (n+ 1J 



where B® = B\ (r — 0). Consequently, 



gS^gg- B % r 2 "- (37) 



We remark that as r — > 0, u dominates over u. From Eq. (jHJ we thus find that the color-flavor locked and unlocked 
condensates coexist in the vortex core and behave as r" near the center. This is similar to the case near H c2 in which 
<f> essentially vanishes. 

We proceed to examine the vortex structure outside the core (r 3> £cfl) by taking into account small deviations 
from the asymptotic behaviors discussed in Sec. lTTTl Su(r) = u(r) — 1, Sv(r) = v(r) — l, and SA%(r) = A^^ — n/gcp-^r. 
The field equations of linear order in these deviations read 

- 2K T (-rf + -4-]su+[a+ (10/3i + 6f3 2 )\n A \ 2 ]5u + S(3i\k a \ 2 5v = 0, (38) 
V dr £ r dr I 



- 2K T (-rt + -4-] Sv + [a + (14/3! + Q(3 2 )\n A \ 2 ]8v + 4/? 1 | Kj4 | 2 <5w = 0, 
\ dr z r dr J 



(39) 



l\l* SA * - as:"- (40) 



It is convenient to diagonalize Eqs. (|38H and l|39() as 

d 2 1 d \ Su + 2Sv 2 5u + 2Sv 



dr^rdrj 3 ~& ^ W 



>CFL 



d 2 1 d \ 2{8u - Sv) 2 2{8u - 5v) 



dr 2 rdrj ^3 e£cFL V3 



(42) 



where e = 3/?cfl//?2- Note that (<5u + 25v)/3 represents a change in magnitude of the CFL condensate from the 
homogeneous one, whereas 2(Su — Sv)/y3 represents the magnitude of the color-flavor unlocked condensate [see Eq. 
©]■ We thus find from Eqs. (|41|l and (|42l) as well as the r —> behaviors discussed above that the magnitude of the 
CFL condensate increases with r in proportion to r n near the vortex center and, above a scale of ^cfl/v^, approaches 
exponentially the homogeneous solution (u = v = 1); the magnitude of the color-flavor unlocked condensate increases 
with r in proportion to r" near the vortex center and, above a scale of ^/e/2£cFL, decreases exponentially to zero. We 
remark that for r ^ Acfl, 5A® v follows the London limit behavior of the supercurrent (|23(l via 5A S V = —n^p^^T \, 
leading to B\ ~ (n</f 8 /27rAQ FL )(7rAcFL/2r) 1 / 2 exp(— r/XcFh)- The vortex structure thus clarified is schematically 
illustrated in Fig. 1. 

The parameter e = 3/?cfl//?2 characterizing the scale of the spatial variation of the color-flavor unlocked condensate 
plays a important role in determining the nature of the interaction between widely separated vortices. For the purpose 
of estimating the intervortex interaction, it is instructive to write down the solutions to Eqs. 14U|) . Q4ip. and l|42|l as 



ln(3CFL,AcFL&A^) - hi 



\CFL 



■0(1), (43) 
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\ 



2(u -v 



FIG. 1: Schematic variation of the mixed field B s z (solid line), the CFL condensate (u + 2v)/3 (dotted line), and the color-flavor 
unlocked condensate 2(u — u)/v3 (dashed line) for a singly quantized vortex (Acfl > v^cfl > £cfl). 



In 



5u + 2Sv 



In 



V2r 



;cfl 



•0(1), 



(44) 



In 



2 (fa - Sv) 
V3 



In 



V2r 



CFL 



0(1), 



(45) 



where K m are the rath modified Bessel functions. Then, following Ref. [l5|], one can estimate the interaction energy 
per unit length of two vortices at large separation a v as 



U(a v ) ~ K t \ka\ 



C 3 K 



CFL 



C 4 ^0 



CcFL 



\/2a v 



CFL 



(46) 



where C3, C4, and C5 are positive dimensionless constants. Since Kq(x) is positive definite and decreases with x, the 
term associated with Acfl in Eq. I|46|l induces repulsion between the vortices, while the terms associated with £cfl 
induce attraction. This reflects the fact that the total magnetic energy of the two vortices becomes larger for smaller 
intervortex separation, while the total loss of the condensation energy arising from the two vortices becomes smaller. 

One of the conditions for type II (I) superconductors is that the intervortex interaction is repulsive (attractive) at 
large separation. In the case of repulsion, a vortex lattice is stable, while in the case of attraction, widely separated 
vortices tend to approach each other more and more closely. From Eq. (|46(l , whose sign is dominated by the exponential 
behavior of the Kq functions [IiiKq(x) = — x + 0(hix) for large x] 7 we find that the intervortex interaction is repulsive 
when kcfl > max(l, y/e)/y/2. We thus see the role played by e in determining the nature of the intervortex interaction. 
Note that e = 4 in weak coupling. 

The type II region is most strongly restricted by the condition that the intervortex interaction is repulsive. In this 
case, the critical value of kcfl is equal to or even larger than l/v2, while the value of kcfl at a s = and that at 
H c = H C 2 are smaller than l/y/2. This is a contrast to the case of ordinary superconductors in which, owing to the 
simple order-parameter structure, all three critical values of the Ginzburg-Landau parameter are degenerate at 1/ y/2. 
We remark that in weak coupling, where kcfl = v/6/7C(3)(67r 2 T c / ff cFLA*) < V2 .10], the CFL condensate near T c 
is deep within the type I region. However, it is still uncertain whether it is type I or type II at low densities where 
T c can be ~ O.l/i and $cfl ~ 1 M- This uncertainty affects the estimates of H c \ and H C 2 since H C 2 = 2kqflH c and 
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V. CONCLUSION 



Within Ginzburg-Landau theory, we have derived the lower and upper critical fields, H c \ and H c i 1 for the CFL 
state, as well as the critical Ginzburg-Landau parameter that distinguishes type I and type II superconductors. In 
doing so, we have examined the structure of a magnetic vortex near the center and outside the core. This vortex is 
associated with the winding of the SU(3) phase in color-flavor space, rather than the U(l) electromagnetic phase, 
and hence allows the presence of a color-flavor unlocked condensate in the core. We have found that the nature of 
the intervortex interaction, controlled by the color-flavor unlocked condensate, provides a stringent constraint on the 
type II regime. 

Application of the present analyses to the magnetic structure and evolution of neutron stars is not straightforward. 
Not only does it require the extrapolation of the weak coupling results near T c to low densities and temperatures, but 
also the clarification of a role of nonze ro q uark masses and neutrality of electric and color charge in modifying the 
CFL condensate at finite temperatures [la. H71 Il8| . Due to uncertainties associated with the extrapolation, it is still 
open whether the CFL condensate, if occurring in stars, would be type I or type II. How nonzero quark masses and 
charge neutrality affect the response to magnetic fields is another open problem. Moreover, as discussed in Ref. [jj, 
to assess the actual situations possible in neutron stars one must take into account the history of the expulsion of the 
magnetic field and the possibility of freezing in of the magnetic field. 
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